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Abstract. The complete near infrared (0.9-2.5 /im) spectra in three different star forming regions (HH24-26, HH72 
and BHR71) are presented and analyzed in the framework of shock excitation models. The spectra are dominated 
by H2 rovibrational emission (vibrational state v < 5, excitation energy T ex < 35 000 K), while emission from 
ionized material, recognizable from [Fe II] and [S II] lines, is significantly fainter. The analysis of the H2 excitation 
diagrams points to the existence of two different excitation regimes: whilst condensations observed only in the 
infrared appear to have temperatures rarely exceeding 3000 K and can be modelled in the framework of steady- 
state C-shock models, the infrared counterparts of Herbig Haro (HH) objects exhibit a temperature stratification 
with components up to more than 5000 K. The H2 emission from representative HH objects (HH26A, HH72A 
and HH320A) has been successfully modelled by planar J-shocks with magnetic precursors, for which the main 
parameters (pre-shock density, speed) are derived. However, these same models are unable to reproduce the 
observed atomic and ionic emission, which probably arises from a distinct and perhaps more embedded region 
with respect to that traced by the H2. Some of the physical parameters of such regions (fractional ionization, 
density) have been estimated in HH72, on the basis of the observed ionic lines. 
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1. Introduction 

The physical effect exerted on the ambient medium by the 
impact of jets from accreting protostars results in the for- 
mation of shock waves. The compressed and heated gas 
radiates away the accumulated thermal energy through 
the emission of ionic, atomic and molecular lines, in rel- 
ative proportions which depend on the structure of the 
shock wave. During the gas cooling, a fundamental role 
is played by molecular hydrogen, whose high abundance 
compensates the low transition rates due to its homonu- 
clear nature. In the near infrared, the association of H2 
emission with knots of shocked gas has been demonstrated 
by mapping at sub-arcsecond scale in the u=l-0 S(l) line 
at 2.122 ^m (e.g. Eisloffel 2000), and the rovibrational 
lines are effective probes of the molecular gas at thousands 
of kelvin, which is missed by optical observations. H2 line 
intensities and intensity ratios are also extensively used to 
distinguish between fluorescence and shock excitation as 
possible mechanisms at the origin of the emission, since 
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strong H2 lines with v>6 in the 1.0-1.4 //m range are ex- 
pected if fluorescence is responsible for the emission (Black 
& van Dishoeck 1987). It has proved more difficult to de- 
rive the nature (C-ontinuous or J-ump type, e.g. Kaufman 
& Neufeld 1996, Hollenbach & McKee 1989, Draine 1980) 
and the physical parameters of the shock waves, such 
as the velocity, the preshock density, the strength of the 
transverse magnetic field and the temperature of the neu- 
tral gas. In the conventional scenario, C-type shocks can 
exist up to velocities typically ss 50 km s _1 , beyond which 
the H 2 molecule is collisionally dissociated and the kinetic 
temperature of the gas rapidly increases, giving rise to 
a discontinuity in the shock parameters. However, more 
recent models (e.g. Smith 1995, Le Bourlot et al. 2002, 
Flower et al. 2003) have shown that H2 dissociation can 
be inhibited over a wider range of shock parameters: from 
slow (v s ^ 25 km s _1 ), partially dissociative J-type shocks, 
up to fast C-type shocks travelling at 70-80 km s _1 . From 
an observational point of view, the presence of different 
shock components (J- or C-type) can be demonstrated by 
the modifications induced in the H 2 excitation diagram: 
the lower post-shock densities attained in C-shocks result 
in a larger departure from LTE conditions than in J-shocks 
of the same speed (Flower et al. 2003). In order to probe 
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such deviations, it is essential to investigate spectroscopi- 
cally the 1.0-1.4 /jm range, since several H 2 lines with dif- 
ferent vibrational quantum number and of relatively high 
excitation energy (> 15 000 K) lie at these wavelengths. 
With the aim of observing such lines, we have undertaken 
a spectroscopic survey (from 0.9 to 2.5 /zm) of a sample 
of Herbig-Haro (HH) objects and H 2 jets, based on the 
observations gathered with the SOFI spectrometer at the 
ESO-NTT. We have already reported the results of this 
survey regarding HH43, HH111, HH240/241 and HH120 
(Giannini ct al. 2002, Nisini et al. 2002, hereafter Paper I 
and II, respectively), showing that quite different excita- 
tion conditions can occur in HH objects: while in HH43 
the bulk of the cooling occurs in H 2 lines, the spectra 
of the other three objects are dominated by iron lines. 
In HH43, the H 2 emission has been successfully fitted by 
means of high velocity, C-type shock models (Flower et al. 
2003), whereas, in the other three cases, the strong ionic 
emission testifies in favour of a dissociative component 
in the shock structure. Here we present the observations 
obtained in three other regions, namely HH24-26, HH72 
and HH320/321 (BHR71). All the spectra exhibit copious 
H 2 rovibrational emission, together, in some cases, with a 
fainter atomic and ionic component. Therefore, they rep- 
resent valid tests of the ability of current shock models to 
predict simultaneously molecular and atomic/ionic emis- 
sion and an opportunity to check the validity of the un- 
derlying assumptions of the models. 

The structure of the paper is the following: we present 
the targetted regions in Sect. 2 and then describe the ob- 
servations and the results obtained (Sect. 3). In Sect. 4, 
we derive the physical parameters of the emitting gas and 
model the observed emission, assuming shock excitation. 
Sect. 5 summarizes our conclusions. 

2. The investigated regions 

2.1. The HH 24-26 region 

Located in the L1630 Orion dark cloud (d « 400 pc, 
Anthony- Twarog 1982), the HH24-26 region is an active 
site of star formation, rich in young stellar objects in 
different evolutionary stages. Two low luminosity, Class 
protostars (Andre et al. 2000), HH24MMS (Chini 
et al. 1993, a 20 oo=05 h 46 m 06.7 s , (5 2 ooo=-00 10'40.8", 
F=205.49,& /J =-14.57, L ho; =5L ) and HH25MMS 
(Bontemps et al.1995, a 200 o=05 h 46 m 06.7 s , 
^ooo=-00°13'24.7", Z J/ =205.53,fe /7 =-14.60, L 6o( =6L Q ) 
drive distinct compact jets, traced by means of 
both the 2.12 /im H 2 line and CO mm transitions 
(Gibb & Davis 1998; Bontemps et al. 1996), while 
the Class I protostar HH26IR (Davis et al. 1997, 
a 20 oo=05 h 46 m 03.9 s , <5 2 ooo=-00°14'52.5", F=205.54, 
b n = — 14.62, Lboi=28.8L ) gives rise to a more extended 
molecular outflow (Gibb & Heaton 1993). Imaging in 
the 2.12 /im H 2 rovibrational line (Davis et al. 1997, see 
Figure 1) reveals a string of shocked H 2 knots along the 
axis of the HH25MMS and HH26IR outflows. Herbig- 



Haro objects 1 (HH24, HH25 and HH26) are associated 
with each of these three sources. Proper motions were 
measured in HH25 and 26 by Chrysostomou et al. (2000), 
who found velocities ranging from less than 70 up to 
200 km s _1 ; the lowest value pertains to HH26A, which 
should be a stationary shock region, and the largest 
one to HH26C in the working surface of a bow shock. 
Far infrared spectra of the region were obtained by 
Benedettini et al. (2000), who interpreted the observed 
emission (from [Oi], CO, H 2 and H 2 lines) as due to a 
mixture of both C- and J-type shocks. Recently, Fabry- 
Perot observations have demonstrated the existence of 
elongated H 2 emission, on scales of a few arcsec, from 
HH26IR and HH25MMS, which appears associated with 
the base of the larger scale jets (Davis et al. 2002). 

2.2. The HH72 region (11660) 

The L1660 region is a bright-rimmed core illuminated 
by the young OB cluster NGC2362 in Vela (Reipurth 
& Graham 1988), at a distance of w 1500 pc (Hilton 
& Lahulla 1995). An E-W, poorly collimated outflow 
(Schwartz, Gee & Haung, 1988) is present in the re- 
gion, possibly driven by the intermediate luminosity 
IRAS source 07180-2356 (Reipurth et al. 1993, a 200 o= 
07 h 20 m 10.3 s , <5 200 o=-24 02'24", F=237.52,& /7 =-4.83, 
L boi w 17OL ). Figure 2 shows the SOFI 2.12 fim im- 
age of the region, where, together with the chain of HH 
objects (HH72 A-C, Reipurth & Graham 1988), several 
pure H 2 features are recognizable (HH72 D-Z). A jet of 
a few arcseconds has been recently observed close by the 
infrared source HH72-IRS (Davis ct al. 2002, sec Fig.2). 

2.3. The BHR71 region 

Figure 3 shows the H 2 2.12/im image by Bourke (2001) 
of the Bok globule BHR71 (Bourke, Hyland & Robinson 
1995), located at w 200 pc (Seidensticker & Schmidt-Kaler 
1989). An embedded binary protostellar system, separated 
by w 3400 AU, has been found in the region (Bourke 2001). 
One of the two sources has been identified as a Class pro- 
tostar (IRS 1: a 2000 = 12 h 01 m 36.6 s , <5 2 ooo=-65 08'48.2", 
Z J/ =297.73, 6 // =-2.78, L 6oi «9L ) by means of 1.3 mm 
observations (Bourke et al. 1997); ISO mid-IR maps led 
to the discovery of the possibly more evolved source IRS 
2 (Bourke 2001, a 200 o = 12 h 01 m 34.0 s , <5 20 oo=-65 08'44.3", 
F=297.72, 6 /7 =-2.78). A recent study by Lemaire et al. 
(2002) suggests the presence of a circumstellar disk around 
this object. A large scale CO outflow, where a shock en- 
hanced chemistry has been revealed (Garay et al. 1998), 
is driven by IRS1 along the North-South direction, while 
IRS 2 gives rise to a more compact outflow with an incli- 
nation m— 36°. Optical [Sll] images revealed HH associ- 
ations (HH320 and HH321) along both the CO outflows 
(Corporon & Reipurth, 1997); these should be at a very 

1 In the following we will refer to HH objects as the optically 
visible condensations. 
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low degree of excitation, as demonstrated by the lack of 
emission in the Ha line. Faint H 2 emission features were 
found by Bourke (2001). 

3. Observations and Results 

The observations were carried out during two runs 
(January 2001 and 2002) with SOFI (Lidman & Cuby 
1999) at ESO-NTT 3.5m telescope (La Silla). Long slit 
spectroscopy was obtained through the blue (0.95-1.64 
/nm) and the red (1.53-2.52 fim) grisms in the low resolu- 
tion mode (A/ A A « 600 with the 1 x 290 arcsec slit). To 
cover most of the H2 emission knots present in the HH24- 
26 region we used five different slit orientations, which are 
depicted in Figure 1. In the HH72 region, we adopted four 
slits to observe the knots both in the blueshifted (knots 
L-B) and redshifted (knots Y and Z) lobes of the outflow 
(Figure 2), while the two jets present in BHR71 were in- 
vestigated with three slits (see Figure 3), covering both 
the Herbig-Haro objects 320/321 and the fainter H2 emis- 
sion features recognizable in the 2.12 fim image by Bourke 
(2001). The observations were performed by nodding and 
jittering the telescope, keeping the target along the slit 
in the usual ABB 'A' mode, with a total integration time 
of 1200 s. The data reduction and calibration were per- 
formed by using the IRAF package. Each observation was 
flat fielded and sky subtracted, while to remove the at- 
mospheric spectral response a couple of telluric O-type 
stars were observed before and after each on-source ob- 
servation. Each target spectrum was then divided by the 
telluric spectrum, corrected by the blackbody function at 
the temperature of the star, and by a few intrinsic hydro- 
gen absorption lines. Wavelength calibration was obtained 
from the spectrum of a Xenon- Argon lamp; the estimated 
wavelength error is smaller than the spectral resolution el- 
ement (2 pixels w 20 A), which corresponds to a velocity 
of about 300 km s _1 . No evidence of shifts larger than this 
limit has been found. Flux calibration (associated uncer- 
tainty 10%) was obtained from the observation of some 
photometric standard stars from the Carter & Meadow 
(1995) catalogue. Ratios between lines lying in the blue 
and red parts of the spectrum, are affected by a 10% un- 
certainty, due to the different instrumental response in the 
two segments. Such error has been estimated by using a 
line (e.g. [Fen] 1.64/xm) present in both parts. As an ex- 
ample of our observations we show in Figs. 4, 5 and 6 
the spectra of HH26A, HH72A and HH320A, which are 
among the knots richest in lines; they are all dominated 
by H2 emission in form of rovibrational lines coming from 
levels with vibrational quantum number up to 5. A rela- 
tively faint atomic and ionic emission is also observed: in 
addition to the bright [C 1] doublet at 0.98 /xm, a few lines 
of [Fe 11] , [S 11] , and [N 1] were detected in the most promi- 
nent knots. Tables la-c, 2a-b and 3a-b (which are available 
only in the electronic version of the paper), list all the de- 
tected lines, along with their spectral identification and 
the vacuum wavelength. Some knots (e.g. HH26A,C) were 
intersected by two slits, but we did not find a significant 



variation in the line ratios along the different directions, 
although the emission is generally brighter along the jet 
axis. To the latter correspond the fluxes reported in the 
Tables. The integrated fluxes were obtained by fitting the 
line shape with a single (or double in case of blending) 
gaussian profile. The reported uncertainty derives from 
the rms of the baseline, multiplied by the width of the 
instrumental profile (30 A). 

The H 2 1-0 S(l) 2.12^m line fluxes observed in HH24-26 
and in HH72 are systematically lower, by a factor between 
1.5 and 10, than the values derived from imaging by Davis 
et al. (1997). The discrepancies are attributable to the spa- 
tial extent of the objects observed, relative to our 1" slit 
width; this explanation is supported by the fact that the 
largest differences occur in those knots, such as HH26D, 
that we observed along their shortest axis (see Figure 1). 
From inspection of Tables la-lc, we see that, along the 
HH26 flow, copious H 2 emission is detected only close to 
the driving source (HH26A and C), while in all the other 
knots faint and weakly excited H 2 lines are revealed. The 
opposite situation exists along the HH72 jet (Tables 2a- 
2b), where the richest H 2 spectra are emitted by the knots 
at the apex of the blueshifted lobe (HH72A and B). Only 
H 2 emission has been detected toward the BHR71 jets 
(Tables 3a-3b). Particularly remarkable is the lack of the 
[Ci] doublet, which, in addition to the present observa- 
tions, has been detected in all the HH objects previously 
investigated (Paper I, II). We also note that, compared 
with the other Herbig-Haro objects in BHR71, HH321B 
has a significantly weaker spectrum, probably because the 
emission peak was badly covered by the slit. 
To complement the SOFI observations, we obtained in 
April 2002 the 3-5 /zm spectra of the HH72 and BHR71 
regions with ISAAC at the VLT (Cuby et al. 2003). In 
this spectral range, several pure rotational lines of H 2 
with high excitation energy (up to « 40 000 K) are lo- 
cated, which can be used to determine the density of the 
molecular gas (Giannini et al. 2001). Unfortunately, the 
poor atmospheric conditions prevented us to obtain high 
quality data at wavelengths larger than 4 /iin, where both 
the thermal emission and the presence of strong water 
absorption bands critically affect the reliability of the ob- 
servations. We used the 1x120 arcsec slit (A/AA w 2000), 
pointing at the Herbig-Haro objects 72 A/B (PA. 84°, i.e. 
slit I in Fig.2) and 320A/B (P.A.320 , i.e. slit I in Fig.3). 
The observational strategy and the procedure adopted 
for the data reduction are similar to those reported by 
Giannini et al.(2001). We detected two H 2 lines: the 1-0 
0(5) at 3.235 fim in HH72B and the 1-0 S(13) at 3.847 
^m in HH72B and HH320A (Tables 2a, 3a). The extracted 
spectra are shown in Figure 7. The 1-00(5) line appears 
narrower of about factor of 2 than the instrumental pro- 
file; however, since its trace in the bidimcnsional image of 
the spectrum is clearly recognizable, we consider the line 
to be a genuine detection. The estimated fluxes were used 
in conjunction with the SOFI data. No intercalibration 
correction was applied to the ISAAC and SOFI lines, but 
the consistency of these two data sets is supported by the 
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Fig. 4. Spectrum between 0.95 and 2.5 ^im of HH26A. 

fact that the column densities of lines coming from the 
same upper level are equal to within the error bars (see 
Sect. 4.2). 

4. Analysis and Discussion 

4.1. Ay determination 

The line lists, reported in Tables 1-3, clearly suggest that 
some differences in the H 2 emission exist throughout all 
the observed knots. In particular, lines coming from the 
highest excited levels (y > 3) are observed mainly in 
condensations visible also in the optical (e.g. HH26A/C, 
HH72A/B, HH320/321), while low excitation conditions 
are recognizable, with few exceptions, in the other knots 
(e.g. HH25, HH72E-Z). Since the most excited lines are lo- 
cated mainly in the blue part of the spectrum (i.e. at wave- 
lengths shorter than 1.4 fim), their absence in the pure in- 
frared knots could be interpreted as being due to a higher 
visual extinction, without necessarily implying lower ex- 



citation conditions. To investigate this crucial point, we 
have firstly estimated the visual extinction in the various 
knots by using pairs of transitions coming from the same 
upper level, for which the difference between the observed 
and the theoretical flux ratio is a function only of Ay (cf. 
Gredel 1994). The accuracy of this method relies strongly 
on the presence of pairs of bright lines at significantly 
different wavelengths which are not affected severely by 
observational problems. In practice, sufficient numbers of 
such pairs of lines are present in the spectra of only the 
brightest, optical knots. For this reason, the values of Ay 
derived in the brightest objects were applied also to the 
nearby knots, where insufficient lines are present to al- 
low a specific determination. Given the large number of 
molecular hydrogen lines, the criteria cited above are sat- 
isfied by several pairs of lines coming from the same upper 
level (e.g. of the kind 2-0 S(i)/2-0 Q(i+2) or 3-1 S(i)/3-2 
S(i)). By using all the available ratios, we were able to 
estimate the visual extinction with an uncertainty of 1- 



T.Giannini et al.: On the excitation of the infrared knots along protostellar jets 



HH72 A 



~i 1 1 r 



to — 

~m m 

— o 

i i 

co cm 




.95 



1.15 



e 

I 

w 

CM 
I 

E 

o 

DO 
CD 



I 

O 



~1 T" 

s 

»3 



~ o* "3 raVcr 
O fc. — o 

I ' — ' n I I 
CO CM 




in I 

o. X 



O CO 

cm 

1 I 

CM f 



1.2 



1.25 



1.3 




- 1 

~n 


i i i i i i i 

* CM CO — 


i i i | i i i i | i i i i | 

S Se S d.o)Qi-Q6 


— r 


i i - 


E m 
— o 
E i 


mm mm 
— o, — o 


So of m m — > 
— CM © — —> 

1 1 1 1 w 








CM — CM — 


CM CO — CM CO 

CM 


1 





2.1 2.2 2.3 

wavelength (/i.m) 



2.4 



Fig. 5. Spectrum between 0.95 and 2.5 fim of HH72A. 

2 mag. Adopting the reddening law of Rickc & Lebofsky 
(1985), we found Ay up to 8 mag in the HH24-26 and 
HH72 flows, while the visual extinction in the BHR71 re- 
gion is small ( Ay < 2 mag), however well compatible with 
the Ay values (< 3 mag) estimated towards the galactic 
plane direction (e.g. Hakkila et al. 1997). The values of Ay 
for individual sources are given in Figs. 8-13. These results 
indicate that the differential extinction only marginally bi- 
ases the observations in all the knots except those having 
Av=8mag (HH25B, HH72J-Z). Thus the presence or ab- 
sence of highly excited H2 lines in the spectra reflects in 
almost all cases the excitation conditions of the molecular 
gas. In this respect, the fact that we detect no lines with 
is>6 indicates that fluorescence mechanisms do not play 
an important role in the excitation (Black & van Dishoeck 
1987). 

An independent determination of Ay could be obtained 
from the [Fell] lines at 1.257 and 1.644 /im. However, 
in our objects, these lines are generally detected with a 



poor signal-to-noise ratio, with the only exception being 
HH72A. For this object, we derive Avpeii]— 4.7±1.5 mag, 
a value higher than the negligible extinction indicated for 
this object by the H 2 lines. The possible implications of 
this result will be discussed below. 

4.2. H2 emission 

In order to derive an estimate of the temperature 
along the flows, we have used the excitation diagrams 
(Boltzmann plots) for the H 2 lines (Figures 8-13), 
plotting the extinction-corrected column densities of the 
rovibrational levels, divided by the statistical weights, 
against their excitation energies. In each plot, lines 
coming from different vibrational states (from v={) to 
v=h) are depicted by different symbols. To minimize the 
uncertainties, we considered only unblended lines with 
a signal-to-noise ratio greater than 3. Any significant 
deviation of the ortho:para ratio from its statistical value 
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Fig. 6. Spectrum between 0.95 and 2.5 /um of HH320A. 
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Fig. 7. H 2 lines observed with ISAAC in HH72B and 
HH320A. 

of 3 would be apparent as a misalignment of the ortho 
and para data points in the rotation diagrams. No such 



deviations are observed: see below and Figs. 8-13. 
If the gas is excited at a single temperature, a straight 
line can be fitted through the data points: this seems 
to be the case only where lines with excitation energies 
lower than typically 15 000 K are detected, i.e. in all the 
molecular hydrogen condensations and few Herbig-Haro 
objects (such as HH321B and HH25A). Here a single 
temperature typically between 2000 and 3000 K can 
account for the observed excitation diagrams. 
Other Herbig-Haro objects (HH26A, HH72A/B, 
HH320A/B and HH321A) and HH26C, which show 
a richer H 2 spectrum, deserve a more realistic approach, 
which consists of relaxing the assumption of a single 
excitation temperature and fitting the various vibrational 
series separately. Although the average temperature 
remains quite similar to that measured in the molecular 
condensations, deviations from the thermal equilibrium 
are evident, with a spread from 2000 K to more than 
5000 K. Similar conditions have been found in two out of 



T.Giannini et al.: On the excitation of the infrared knots along protostellar jets 



HH 26 



HH72: HH objects 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 
T v =3000±160 K; T v _ v> =3970±210 K 

1 =3420±110 K 

,4^2 mag 

HH26A 

I I I I | I I I I | I I I I | I I I I | I I I I | 
5 T V =2590±110 K; T„=4510±640 K; T =5170±630 K 

3350±150 K . 
^52 mag 





| 34 



M 32 
Z 

30 

28 
36 



no 32 
Z 

a 30 



1 1 1 1 1 1 1 
: HH72A 


l v — cD4U± loU 1\, 1 ¥ _ T 


i i i i | i 
— iRpri-i-TRn v 

-oooU±JOU J\ 

-29foO±120 K _ 
A,^2 mag 


I I I I | I I 


I I I I I I 

T v =2290±160 K; T v __ v 


~-l 1 1 1 | 1 
=4960±690 K ! 




T 

me 


=3150±220 K : 






A,S2 mag ~. 


: HH72B 








i l i i i i l i i i i 


i ~r-i-. ,1," 



15 20 
E v , (10= K) 



Fig. 8. Excitation diagram for HH26 A and B. Lines com- 
ing from different vibrational levels are indicated with dif- 
ferent symbols: filled circles, open circles, filled triangles, 
open triangles, filled squares and open squares indicate 
lines coming from levels from to 5, respectively. The 
straight line(s) represents the best fit(s) through the data, 
at the temperature(s) reported in top right angle. The de- 
rived value of the extinction is also indicated. 

HH24-25-26 



Fig. 10. As in Figure 8 for HH72 A and B. 
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Fig. 9. As in Figure 8 for the pure H2 knots in the HH24- 
26 region. The two crosses in the HH25C panel indicate 
the positions which would correspond to the 3-1 S(5) and 
4-2 S(5) lines by assuming the same ratios with the 1-0 
S(l) line as in the HH26C case (see text). 
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Fig. 11. As in Figure 8 for the pure H 2 knots in the HH72 
region. 

the four HH objects we investigated in previous studies, 
namely HH43 (Paper I) and HH120 (Paper II). 
Since the lack of a temperature stratification is deduced 
only when the highly excited lines (with v > 3) remain 
undetected, it is crucial to understand whether, in these 
cases, the observations reflects an intrinsic property of 
the gas or is attributable to the instrumental detection 
limit. 

In order to investigate this point, we have estimated 
the fluxes expected for the v > 3 lines, both from a single 
temperature gas and in the case of a temperature strat- 
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Fig. 12. As in Figure 8 for HH 320A/B and HH321A/B. 
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Fig. 13. As in 

BHR71 region. 



Figure 8 for the pure H 2 knots in the 



ification. As a test case, we considered the spectrum of 
HH25C, which is a bright knot with a negligible extinc- 
tion, observed in only j/=l,2 lines, which can be fitted by 
a single temperature w 2300 K (Figure 9, middle right 
panel). By extrapolating the straight line up to energies 
~ 25 000 K, we can evaluate the fluxes of the lines 
that would be expected in thermal equilibrium: we find 
that all the i^=3,4 lines are definitely below the detection 
limit. We have estimated the fluxes of the same lines, in 
case of a temperature stratification, by assuming that their 
ratios with the 1-0 S(l) line are the same as those mea- 



Fig. 14. Top panel: [Fe n] 1.644 /xm/[S n] 1.03 /im as func- 
tion of the temperature, for an electron density of 5 x 10 4 
cm~ 3 . The curves refer to different assumptions on the 
iron abundance. The observed ratio in HH72A is drawn 
as a hatched area. Bottom panel: as in the top panel for 
the ratio [Fen] 1.644 /xm/[Ni] 1.04 /im. The parameter g 
is the ratio between the iron abundance and the fraction 
of neutral nitrogen. 



sured in HH26C (a knot with brightness and extinction 
comparable to HH25C, but showing a spectrum richer of 
z/=3,4 lines). In Figure 9, the two crosses indicate the cor- 
responding column densities of the 3-1 S(5) and 4-2 S(5) 
lines: both these points lie above the 3tr detection limit. 
Thus we conclude that the absence of these lines in the 
spectrum of HH25C points to excitation conditions differ- 
ent from HH26C. 

This analysis has been repeated for other bright knots 
(such as HH72D and knots 4 and 9 in BHR71), yielding 
the same conclusion as in the case of HH25C. Thus tem- 
peratures rarely exceeding 3000K appear to be a proper 
feature of the pure infrared condensations, while in HH 
objects significant temperature gradients occur. This ten- 
tative conclusion needs to be confirmed by enlarging the 
observational database. 



4.3. Atomic and Ionic emission 

As we have already pointed out, faint atomic and ionic 
lines are detected in only four of the objects investigated 
(HH24A, HH26A and HH72A/B). We first attempt to an- 
alyze such emission without making assumptions regard- 
ing the underlying excitation mechanism, using simple 
NLTE models of the equilibrium level populations. We 
then present the results obtained by modelling the same 
emission, assuming shock excitation (Sect. 4.4.6). 



T.Giannini et al.: On the excitation of the infrared knots along protostellar jets 



9 



Following the procedure described in Paper II, 
we determine the electron density from the [Fen] 
1.644 /tm/1.600 /im ratio, obtaining n c of approximately 
5 x 10 4 cm~ 3 for HH24A and HH72A, a value which can 
be considered also as a 3tr upper limit for HH26A. Since 
all the transitions of [Fe n] detected come from levels with 
quite similar excitation energies, their relative intensities 
are insensitive to the temperature. However, some con- 
straints on the conditions in the ionized gas can be derived 
by using lines from different species, such as the [S n] and 
the [Ni] multiplcts at 1.03 /jm and 1.04 /im. In practice, 
this can be done for HH72A, the only object where these 
lines are all detected with a sufficiently high signal-to-noise 
ratio. 

From the [Fen] 1.644 /xm/[Sn] 1.03 /im ratio we derive 
the fraction of iron in the gas phase, a parameter regu- 
lated in HH objects by the shock efficiency in destroy- 
ing the cores of the dust grains in which most of the 
iron is to be found (cf. May et al. 2000). To this end, 
we plot, in the upper panel of Figure 14, the expected 
[Fen] 1.644 /xm/[Sn] 1.03 /zm ratio as a function of the 
temperature, having fixed the electron density at the de- 
rived value of 5 x 10 4 cm~ 3 . The hatched area refers to 
the HH72A observations, corrected by Av=4.7 mag. We 
assume that all the sulphur is in the gas phase and that 
both species are singly ionized, as suggested by the ab- 
sence from the spectra of transitions from any other ion- 
ization stage. Then, the plotted ratio is a function only of 
the temperature and the fraction /of iron in the gas phase. 
The observational data agree with / values less than 0.4 
for all the considered temperatures (2000 < T < 15 000 
K); in particular, if T=10000 K, we find 0.20 < / < 0.35. 
For comparison, the same analysis applied to HH240A, the 
object in the survey with the richest iron spectrum, gives 
0.6 < / < 1, in substantial agreement with the abundance 
estimated in Paper II by comparing with the hydrogen re- 
combination lines. 

Once the iron abundance has been evaluated, it is possible 
to use the [Fen] 1.644 /xm/[Ni] 1.04 /im ratio to estimate 
the hydrogen ionization fraction (x c ) and hence the total 
density of the ionized medium (nn = n e /x e ). The depen- 
dence of the above ratio on x c assumes that nitrogen is at 
least partially ionized, as is observationally testified by the 
frequent detection in HH objects of the [Nil] optical lines 
at 6548 and 6583 A (e.g. Bally & Reipurth 2001). The 
considered ratio is plotted in the lower panel of Figure 14 
for different values of the parameter g, defined as the ra- 
tio between the iron abundance and the fraction of neutral 
nitrogen. 

If we take T=10 000 K, we derive 0.05 < g < 0.3, which, 
for the estimated range of/, corresponds to N°/N > 0.67. 
Following the procedure described by Bacciotti & Eisloffcl 
(1999) for deriving the nitrogen ionization equilibrium, 
N°/N > 0.67 corresponds to x c < 0.25, and to n H > 
2.0 x 10 5 cm~ 3 . The same procedure applied to HH240A 
gives 0.85 < x c < 0.95, in agreement with the quite high 
excitation/ionization conditions derived for this object. 



4.4. Comparison with predictions of shock models 

In this section, we compare the emission from the observed 
regions with the predictions of shock models. As repre- 
sentative cases, we have chosen the knots richest in lines: 
HH26A, HH72A and HH320A. We also model HH25C, 
which is the 'test case', considered in Section 4.2 above, 
of an object which is observed in only v=l,2 lines. 

The main features and the methods used in the shock 
model computation are described by Le Bourlot et al. 
(2002) and Flower et al. (2003). Here, we summarize only 
those aspects most relevant to the present study. The code, 
MHD_VODE, simulates one-dimensional, planar, multi- 
fluid shock waves. It solves the magnetohydrodynamical 
equations for the three fluids (neutral, positively and neg- 
atively charged) in parallel with an extensive chemical net- 
work, which links 132 species by 917 reactions. 

The manner in which the emission by molecular hy- 
drogen is calculated is also described by Le Bourlot et 
al. (2002) and Flower et al. (2003). Briefly, the processes 
which populate and depopulate the rovibrational levels 
of H 2 are: collisional excitation and de-excitation; spon- 
taneous radiative decay; collisional dissociation and ion- 
ization; and reformation on grains, which occurs in the 
wake of the shock wave. The equations for the H2 level 
populations are solved in parallel with the chemical and 
dynamical conservation equations. This approach is es- 
sential to ensure the accuracy of the computed H 2 column 
densities, because the level populations do not respond in- 
stantaneously to changes in the physical state of the gas. 
Radiative pumping of H 2 was not taken into account and 
is considered to be unimportant, compared with collisional 
excitation, in the context of the models presented below 
and for reasons that will be given there. 

In the presence of a transverse magnetic field, a 'jump' 
(J-) type shock wave will tend to evolve into a 'continu- 
ous' (C-) type shock wave in which, prior to the attain- 
ment of steady state, a J-type shock remains embedded 
(Pineau des Forets, Flower & Chieze, 1997; Smith & Mac 
Low, 1997; Chieze, Pineau des Forets & Flower, 1998). 
MHD_VODE can simulate not only steady state J- and 
C-type shock waves, but also J-type shocks with a mag- 
netic precursor, as well as the equilibrium conditions of 
the preshock gas (see Flower et al. 2003). 

In addition to being observable in rovibrational transi- 
tions of H2, HH outflows can be traced through the emis- 
sion of CO and other molecules and of atoms and ions. In 
the context of the model, as the temperature and density 
of the gas increase, H 2 is dissociated and H begins to be 
ionized. As a consequence, the contribution of radiative 
cooling by H 2 decreases and that of other molecular and 
atomic species, such as H 2 0, CO, O, C, C + and Fe + , 
assumes greater importance. MHD_VODE incorporates 
the cooling due to rovibrational transitions of molecules 
and electronic transitions of atoms, specifically C, N, O, 
Si, S, and of their ions, as well as Fe + . The strongest 
ionic/atomic lines observed are those of [Fen] and [Ci]; 
their predicted intensities are affected strongly by pro- 
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cesses within the shock wave. In a C-type shock wave or 
a magnetic precursor to a J-type shock, the erosion of 
grains can occur owing to the streaming of charged grains 
through the neutral gas. At sufficiently high shock speeds, 
the (refractory) grain cores can be partially eroded, re- 
leasing elements such as Fe and C into the gas phase: 
see May et al. (2000). Neutral iron is ionized rapidly in 
the gas phase, predominantly through charge transfer re- 
actions with ions, such as H3 + and H30 + , which have 
larger ionization potentials. Fe + may then be excited col- 
lisionally and emit the [Fe n] forbidden line spectrum. The 
predictions of the model with regard to these two species 
are considered in Sect. 4.4.6 below, and information on 
the atomic data used in the model can be found in the 
Appendix. 

4.4.1. Model Parameters 



Table 4. Parameters of shock excitation models 



Object 


v s 




Age 


Radiation field 




km s _1 


cm" 3 


yr 




HH72A 


31 


10 4 


185 


Yes 


HH26A 


52 


10 4 


70 


Yes 


HH320A 


41 


10 4 


475 


No 


HH25C 1 


45 


10 4 


> 260 


Yes 


HH25C 2 


27 


10 5 


> 680 


Yes 



Notes: a n H = n(H) + 2n(H 2 ) 
\ 2 Models 1 and 2 of HH25C 



It was found that, for HH72A, HH26A and HH320A, 
neither steady state C-type nor steady state J-type shocks 
(which are obtained when the magnetic field is sufficiently 
weak) could satisfactorily reproduce the observations. The 
former underestimated the column densities of the levels 
of high excitation, whereas the latter underestimated the 
column densities of the levels of low excitation. However, 
non-equilibrium models of J-shocks with magnetic precur- 
sors provided good fits to the H2 observations of these 
three HH objects. On the other hand, the excitation di- 
agram of HH25C could be reproduced by a steady state 
C-type shock model. In this case, the introduction of a 
J-component resulted in the populations of levels of high 
excitation being overestimated. 

The optimal values of the shock speed and preshock 
gas density (in the range from 10 3 to 10 5 cm~ 3 ) were 
determined by comparing the predicted H2 excitation 
diagram with the observations. The relation B(pG) = 
[nnfcm -3 )] - 5 was adopted as defining the initial value 
of the magnetic induction. Two additional parameters are 
the age of the shock wave, which determines the extent of 
the magnetic precursor, and the strength of the ultraviolet 
radiation field. The model parameters that produce opti- 
mal fits to the excitation diagrams for each object are pre- 
sented in Table 4. The entry in the final column indicates 



whether the effects of an ultraviolet radiation field (via 
photoionization processes) were included in the model. 
HH72A, HH26A and HH25C arc located in the vicinity of 
sources of ultraviolet radiation, and, for these objects, we 
included photoionization processes with rates correspond- 
ing to the mean interstellar radiation field in the solar 
neighbourhood (Draine, 1978). On the other hand, there 
appear to be no sources of ultraviolet radiation nearby 
HH320A, and, in this case, photoionization processes were 
neglected. It is perhaps significant that the inclusion of a 
radiation field improved the fits to the empirical excita- 
tion diagrams for HH72A, HH26A and HH25C but made 
worse the fit for HH320A. 

The rate of photodissociation of H 2 by the mean in- 
terstellar radiation field is of the order of 10~ n s _1 (cf. 
Abgrall et al. 1992). The same process (of photon absorp- 
tion) which gives rise to photodissociation also gives rise to 
H2 fluorescence through radiative cascade from the excited 
(optically pumped) electronic states to the rovibrational 
levels of the electronic ground state. Optical pumping fol- 
lowed by fluorescence is an order of magnitude more prob- 
able than optical pumping followed by dissociation of H2. 
Thus, the total rate of the fluorescence cascade is of the 
order of 10~ 10 s _1 . On the other hand, any given rovibra- 
tional level of the electronic ground state receives only a 
fraction of the cascades. We take 10~ n s _1 as an order of 
magnitude estimate of the rate of population of any given 
excited rovibrational state through the optical pumping 
by the unshielded mean interstellar radiation field. We 
shall now compare this rate with that of collisional exci- 
tation. 

At kinetic temperatures of a few thousand kelvin, at- 
tained in the shock wave, the rate coefficients for rovi- 
brational excitation of H2 by H are typically of the order 
of 10 -12 cm~ 3 s _1 and exceed those for excitation by He 
and H2 by more than a factor of 10 (cf. Le Bourlot et 
al. 1999). Furthermore, the fractional abundance of H in 
the shock wave becomes comparable with that of He, i.e. 
of the order of 0.1. In the models considered below, the 
preshock gas density nu > 10 cm~ 3 . Thus, the rate of 
rovibrational excitation of H2 by H is of the order of 10~ 9 
s _1 or greater. We conclude that the inclusion of optical 
pumping would not have a significant effect on computed 
column densities of the excited rovibrational levels of H 2 
and hence on their comparison with the observations. 

Low preshock gas densities result in low rates of col- 
lisional excitation of H2 molecules, smaller excited state 
column densities and greater departures of the level pop- 
ulations from LTE. For a given density, the shock speed 
is the main factor determining the column densities: an 
increase in the shock speed results in larger column den- 
sities. Indeed, the model is so sensitive to the shock speed 
that a change of the order of 1 km s _1 in f s can make a 
significant difference to the excitation diagram. The age 
of the outflow affects differentially the population den- 
sities of the more highly excited levels: the younger the 
object, the closer is the shock wave to being J-type and 
the higher are the column densities of the more highly ex- 
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cited rovibrational levels. Furthermore, the populations of 
these levels tend to be closer to LTE. The presence of a 
radiation field leads to an increase in the degree of ioniza- 
tion of the gas and hence in the strength of the coupling 
between the neutral and the charged fluids; there results 
a narrower shock wave. 

The observations comprise emission lines from H2 lev- 
els of both ortho (J odd) and para (J even) symmetry and 
hence reflect the variation of the ortho:para ratio through 
the shock wave. In the hot gas produced by the shock 
wave, para-H 2 is converted to ortho-H 2 by nuclear spin 
changing reactions with H. As shown by Wilgcnbus et 
al. (2000), the local ortho:para ratio increases from a low 
value in the cold preshock gas (assuming that it has be- 
come thermalized through proton transfer reactions with 
H + and H^~) to a value of 3 immediately behind the shock 
wave. The ratio of the total column densities in the ortho 
and para levels, N(ortho)/N(para), lags behind the local 
value of the ratio of the number densities of the ortho 
and para levels. However, our observations do not yield 
N(ortho)/N(para), as the contributions from the levels of 
the vibrational ground state, whose populations become 
significant in the cooling flow, are not included. Rather, 
the observations (of excited vibrational levels) reflect the 
value of the ortho:para ratio in the region of the hot gas in 
which the corresponding levels are preferentially excited. 
Thus, we consider that any attempt to derive a unique 'or- 
tho:para ratio' from our observations would be misguided. 
On the other hand, any systematic departure of the or- 
tho:para ratio from its value in LTE would manifest itself 
as a downwards displacement of the ortho levels in the ex- 
citation diagram (cf. Wilgenbus et al. 2000); no such shift 
is discernible in the comparison with the models presented 
below. 

4.4.2. HH72A 
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Fig. 15. Velocity and temperature profiles of the model of 
HH72A. The velocity is expressed in the shock frame and 
the preshock gas is to the left. The J-shock discontinuity 
is evident in both the temperature and velocity profiles. 

The velocity and temperature profiles of the model for 
HH72A are plotted against flow time, t n = f —dz n , in 
Figure 15: the time at which the J-type discontinuity oc- 



curs may be identified with the age of the shock wave. 
The models of HH26A and HH320A exhibit qualitatively 
similar velocity and temperature profiles, which will not 
be shown. 

In Figure 16 (top panel) we show the excitation di- 
agrams derived from the model (open symbols) and the 
observations (corresponding filled symbols) . The observa- 
tions extend up to approximately 2 x 10 4 K above ground 
and include ^=1, 2 and 3 vibrational levels. At successive 
vibrational thresholds, the observational points do not de- 
part notably from the median curve, indicating that there 
is a high degree of thermalization at the local temperature. 
Models with initial densities less than 10 4 cm~ 3 yield lower 
column densities which diverge more evidently from the 
median curve at each vibrational threshold. Conversely, 
models with initial densities higher than 10 4 cm~ 3 pro- 
duce level populations which are closer to LTE, but the 
curvature of the Boltzmann plot is not consistent with 
that observed. Similar conclusions were reached from the 
analysis of the spectra of both HH26A and HH320A. 

The characteristic length of the shock wave is deter- 
mined by the distance over which the ion and neutral ve- 
locities differ. For HH72A, this length is approximately 2 
x 10 16 cm, which may be compared with the observed 
dimension of the source, 3.5 x 10 16 cm (at a distance of 
1500 pc and for an angular size of 1.5 arcsec, as measured 
on the image taken at the telescope). The same compar- 
isons for HH26A and HH320A produced similar levels of 
agreement. 

4.4.3. HH26A 

The observed points in the H2 excitation diagram extend, 
in this case, to almost 3 x 10 4 K and include levels of the 
first four excited vibrational states: see Figure 16, middle 
panel. The column densities are larger than observed in 
HH72A, implying a higher shock speed, and show a similar 
degree of departure from LTE. The higher shock speed 
used for HH26A results in a higher rate of H 2 cooling and 
so the shock is narrower than HH72A; it is for this reason 
that the model of HH26A has a smaller age than HH72A. 

The model of HH26A yields column densities of levels 
with v = 1 which exceed the observations beyond approx- 
imately 10 4 K (cf. Figure 16); this is symptomatic of a 
more general tendency of the models to overestimate the 
populations of the levels in the v = 1 manifold. Such a 
tendency may point to an overestimation of the rates of 
vibrationally inelastic v = — > v = 1 collisions. The model 
of HH26A has the highest shock speed of the three ob- 
jects that we consider and hence has the largest fractional 
abundance of H, owing to collisional dissociation of H 2 . As 
the rate coefficients for vibrational excitation of H2 by H 
are much larger than for excitation by He or H 2 , a larger 
H abundance exacerbates the effect of any systematic er- 
ror in the H-H 2 collisional rate coefficients. However, we 
stress that there is no independent evidence to confirm 
the existence of such a systematic error. 
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Fig. 16. Models predictions for the excitation diagrams of HH72A, HH26A and HH320A. Different symbols distinguish 
transitions coming from different manifolds. 



4.4.4. HH320A 

Levels of the first six vibrational states, with excitation 
energies extending to over 3 x 10 4 K, are observed: see 
Figure 16, bottom panel. The column densities observed 
for two levels with v = 3 (E = 20 857 K and E = 23 070 K) 
fall significantly below the mean curve in the excitation 
diagram. The level at E = 20 857 K has been observed 
twice, in two different transitions: one observation falls on 
the median curve, while the other, as already stated, falls 
below it and out of line with the remaining v = 3 levels. We 
believe that this discrepancy arises from systematic errors 
(related to atmospheric absorption) in the observations 
and that the misalignment of the observational point at 
E = 23 070 K level is attributable to a similar effect. The 
absence of a radiation field in this model results in a wider 
shock wave than is the case for the other objects and hence 
to a greater shock age. 

No atomic or ionic emission is observed towards the 
BHR71 region and the absence of both [Fen] and [Ci] 
lines, which were observed towards all other outflows con- 
sidered in this paper and in papers I and II, is particu- 
larly surprising; here we discuss why these lines are not 
observed. 



Perhaps the most extreme possibility is that there is 
very little Fe and C present in the gas phase of the BHR71 
region because they are mainly in solid (and, in the case 
of C, in molecular) form, which would imply that the dust 
grains have undergone little erosion. Whilst such a situ- 
ation cannot be discounted by our observations, we con- 
sider that it is unlikely because, according to our mod- 
elling results, the gas has been subjected to a strong 
shock. Emission from [Ci] and [Fen] would not be ob- 
served if C and Fe + were ionized by collisions, photons 
or charge transfer reactions; but it should be noted that 
we observe no lines of Fe 2+ in this region. The energies 
required for ionization of C and Fe + (11.26 eV and 16.16 
eV, respectively) are such that, if collisional or photoion- 
ization is invoked, we should expect emission from other 
ions (for example, [S II] and [N II]) that are not observed. 
Charge transfer ionization of C to C + is ineffective and 
Fe 2+ rapidily recombines via charge transfer recombina- 
tion reactions with H (Kingdon and Ferland, 1996). For 
these reasons, we consider that ionization is unlikely to 
be responsible for the absence of [C i] and [Fe n] emission. 
Collisional deexcitation of the lines could be responsible 
for their relative weakness. The rate of electron collisional 
deexcitation of the [Fen] and [Ci] lines is of order 10 -7 
and 10 -9 cm 3 s _1 , respectively: by comparison, the radia- 
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tive decay rates are of order 10~ 3 s _1 for [Fen] and 1CT 4 
s _1 for [Ci]. Therefore, in order for collisional deexcita- 
tion to dominate radiative decay, n e must be greater than 
10 or 10 5 cm -3 , respectively, which would imply that hy- 
drogen is fully ionized. An alternative and more probable 
explanation is that the atomic and ionic emission arises in 
the apex of a bow shock or a reverse shock in the jet. In 
either case, the region of atomic and ionic emission would 
be small, and it is possible that the slits did not encompass 
fully the emitting areas. 

4.4.5. HH25C 

As noted in Sect. 4.2, HH25C is an object which is ob- 
served in 1^=1,2 lines, at excitation energies extending up 
to « 1.5 x 10 4 K, beyond which we have only upper lim- 
its on the line intensities and the column densities of the 
corresponding levels. 

We present two models of HH25C, whose parameters 
are given in Table 4. Both the models are of C-type shocks 
in steady state, and the flow time through the shock wave 
provides only a lower limit to its age. The models are 
compared with the observations of HH25C in Figure 17. 
Regarding the levels which are definitely observed, with 
excitation energies less than about 1.5 x 10 4 K, the model 
with the lower density (and higher shock speed) provides 
perhaps a better fit: the column densities of the v=\ levels 
are overestimated and those of the v=2 levels are under- 
estimated to a greater extent by the higher density (lower 
velocity) model 2. On the other hand, the predictions of 
the higher density model are more consistent with the ob- 
servational upper limits at excitation energies above 1.5 
x 10 4 K. Straight line fits to the theoretical data yield 
excitation temperatures of 2770 K for model 1 and 2130 
K for model 2; the latter result is more consistent with the 
value of w 2300 K obtained by fitting to the observational 
data. 

4.4.6. [Fen] and [Ci] emission lines 

Emission lines of [Fen] and [Ci] were also observed to- 
wards HH72A and HH26A. For both objects, the inten- 
sities predicted by the models (which successfully repro- 
duce the H 2 excitation diagrams) are more than 2 orders 
of magnitude (in the case of [Fc n] ) and 5 orders of mag- 
nitude (in the case of [C i] ) smaller than are observed. 

An important difference between the shock wave model 
and the NLTE model used in Sect. 4.3 is that the shock 
model predicts a low ionization fraction. For example, the 
ionization fraction within the shock wave model of HH72A 
does not exceed about 2 x 10~ 5 . This degree of ionization 
is insufficient for electron collisions to play a significant 
role. The excitation of Fc+ and of C occurs principally in 
collisions with neutral perturbers, which are less effective 
than those with electrons. For the rates of collisional exci- 
tation of [Fe n] by electrons to be comparable with the rate 
of excitation by neutrals, the ratio x c — n c /«H ~ 10~ 3 . 



In order to enhance the computed [Fell] line intensities 
to values comparable with those observed requires x « 
0.1; this may be compared with the upper limit x c < 0.25 
derived in Sect. 4.3. 

The fact that the models underestimate the [Fe II] and 
[C i] line intensities suggests that these forbidden lines are 
emitted from regions where the degrees of dissociation and 
ionization are much higher, which do not contribute signif- 
icantly to the H2 line intensities. Such regions could be the 
apex of a bow shock (with the H2 emission being produced 
in the wings), a reverse shock in the jet, or photoionized 
gas. However, the last of these three possibilities would im- 
ply local sources of ultraviolet radiation, sufficiently hard 
and intense to photoionize hydrogen; we know of no ob- 
servational evidence for the existence of such sources of 
radiation. Accordingly, we lean towards an explanation in 
terms of the collisional ionization of hydrogen, in high ve- 
locity J-typc shocks. Indirect evidence that the [Fe 11] and 
H2 emission lines do, indeed, arise in distinct regions is 
afforded by the difference in the values of the extinction 
which are derived from the respective line intensities (at 
least in the case of HH72A) and from the high values of 
the electron density derived from the [Fell] lines. In fu- 
ture studies, we intend to address the issue of the origin 
the atomic and ionic forbidden lines in a more quantitative 
manner, within the framework of a shock model. 

5. Conclusions 

We have presented the 0.95-2.5 /im spectra of three proto- 
stellar jets (HH24-26, HH72, HH320-321). The analysis of 
the observed features has led to the following conclusions: 

- Strong H2 emission lines with excitation energies up 
to 35 000 K are detected throughout all the Herbig- 
Haro objects present in the three regions, while in- 
frared knots show lines with excitation energies rarely 
exceeding 15 000 K. This difference in the line emission 
reflects different excitation regimes: while the conden- 
sations observed only in the infrared are excited at 
a single temperature of « 2000-3000 K, Herbig-Haro 
objects have a temperature stratification, with compo- 
nents up to more than 5000 K, which tend towards the 
values (of the order of 10 4 K) determined from optical 
transitions. We note that, in order to trace effectively 
the highest temperature components, it is essential to 
obtain observations in the 1.0-1.4 /zm range, where the 
lines with high vibrational quantum numbers are lo- 
cated. 

- Atomic and ionic emission in form of [Fe 11] , [S 11] and 
[Ni] transitions is detected in only a few HH objects. 
The ratios of lines of different species have been used to 
derive some of the physical parameters of the ionized 
gas (electron density, iron abundance and ionization 
fraction), which all indicate that the shocks present 
in these objects have lower excitation/ionization con- 
ditions than in cases where stronger [Fen] emission 
occurs. 
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Fig. 17. Excitation diagram of HH25C: the observational data are compared with predictions from two different 
steady-state C-shock models (1 and 2) whose parameters are listed in Table 4. The arrows indicate 3cr upper limits 
on the column densities of the corresponding transitions. 



- The observed H2 emission in the HH objects can be 
reproduced by models of J-type shocks with magnetic 
precursors and ages of typically a few hundred years. 
The shock speeds are in the range from about 30 to 
about 50 km s _1 , and the preshock gas density is of 
the order of 10 4 cm -3 . In the case of the pure H2 knot 
HH25C, it is possible to fit marginally the observa- 
tional data by means of a steady state (C-type) shock 
model with a higher preshock density (10 5 cm -3 ). 

- The same planar model that fits the H2 emission sys- 
tematically underestimates both the [Fell] and [Cl] 
lines fluxes observed in HH72A and HH26A. We be- 
lieve that either a bow shock (with atomic and ionic 
emission originating at the apex of the shock and with 
H2 arising in the wings) or a reverse shock in denser 
gas at the centre of the outflow is responsible for the 
ionic and atomic emission. Supporting evidence comes 
from the high electron density (« 5 x 10 4 cm~ 3 ) and 
visual extinction derived from the [Fe 11] lines. 
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6. Appendix: Fe + and C° models. 

In the model of Fe + , we considered transitions among the 
19 energy levels which arise from the a 6 D, a 4 F, a 4 D, a 4 P 
and b 4 P terms. All of the observed infrared transitions 
are emitted from a 4 D; the strongest lines are at 1.644 /zm 
(a 4 D 2 - a 4 Fs) and 1.257 ^m (a 4 D 2 - a 6 D;i). In order 

v 2 2 ' ^2 2 ' 

to ensure that the populations of the levels of the a 4 D 
term, which lies about 1 eV (w 10 4 K) above ground, are 
converged, we included cascades from the a 4 P and b 4 P 



terms; the b 4 P triplet lies approximately 3 eV («3 x 10 4 
K) above ground. 

The rates of collisional and radiative transitions de- 
termined the Fe + level populations at each point in the 
model. The spontaneous radiative transition probabilities 
computed by Quinet, Le Dourneuf & Zeippen (1996) were 
adopted. Collisions with both neutral species and elec- 
trons were included; for the latter, we used the collision 
strengths of Zhang & Pradhan (1995). Regarding the neu- 
tral perturbcrs, we included collisions with the most abun- 
dant species, H2, H and He. However, we are not aware 
of any quantum mechanical calculations of cross sections 
for the electronic excitation of Fe + by these perturbers. 
Accordingly, we estimated the rate coefficients using a 
simple classical approximation, which derives from the 
Langevin 'orbiting' model, adopting the polarizability and 
reduced mass appropriate to each perturber. 

Assuming LS coupling (cf. Nussbaumer & Storey 
1980), the conservation of the total electron spin implies 
that the spin state of Fe + can be changed only by ex- 
change of a bound and an incident electron with oppo- 
sitely directed spins [if (weak) magnetic interactions are 
neglected]. Electron collisional excitation of the a 4 D term 
from the a 6 D ground term proceeds via this electron ex- 
change mechanism. In collisions between Fe + and H, a 
similar process of exchange can, in principle, take place, 
between one of the bound electrons of the ion and the 
bound electron of the hydrogen atom. In the case of col- 
lisions with He, on the other hand, the ground state is 
a singlet, and exchange with an electron of opposite spin 
implies a transition to a triplet state. As the lowest energy 
triplet state of He lies almost 20 eV above ground, such 
events will be extremely rare under the physical conditions 
that we are considering. Finally, ground state H2 is also 
a singlet, and the lowest energy triplet state is repulsive. 
Hence, transitions to the triplet state lead to dissociation 
of the molecule and require an energy input of more than 
4.5 eV. We conclude that transitions involving a change 
of spin, including excitation of a D from a 6 D, are induced 
at a significant rate only by collisions with electrons and 
H atoms. 

Collisional and radiative transitions among the five en- 
ergy levels of neutral C which arise from the 3 P, 1 D and 
1 S terms were incorporated in the model. Following the 
same reasoning as above, we should include all levels of 
C° which contribute significantly, by direct collisional ex- 
citation or cascade, to the rates of population of the levels 
from which transitions are observed. The uppermost 1 S 
term has an energy of 3.1 x 10 4 K, relative to ground. 
The 0.983 ^m and 0.985 /im lines originate in the 1 D 
term, which has an energy of 1.6 x 10 4 K. Population 
transfer in collisions with H (Launay & Roueff 1977), H2 
(Schroder et al. 1991), He (Staemmler & Flower 1991), 
and H+ (Roueff & Le Bourlot 1990) was taken into ac- 
count. The data pertaining to electron collisions and also 
the radiative transition probabilities were taken from the 
compilation of Mendoza (1983). 
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